Infrared dielectric properties of low-stress silicon nitride 
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Silicon nitride thin films play an important role in the realization of sensors, filters, and high-performance 
circuits. Estimates of the dielectric function in the far- and mid-infrared regime are derived from the observed 
transmittance spectra for a commonly employed low-stress silicon nitride formulation. The experimental, 
modeling, and numerical methods used to extract the dielectric parameters with an accuracy of approximately 
4% are presented. © 2012 Optical Society of America 
OCIS codes: 310.3840, 310.6188, 310.6860. 



The physical properties of silicon nitride thin films, 
namely low tensile stress, low thermal/electrical con- 
ductance, and its overall compatibility with other com- 
mon materials, have facilitated its use in the micro- 
fabrication of structures requiring mechanical support, 
thermal isolation, and low-loss microwave signal propa- 
gation (e.g., [l}|4j)- Silicon nitride films are amorphous, 
highly absorbing in the mid-infrared (5 , and their gen- 
eral properties are functions of composition Ryul- Here, 
the optical properties are studied in detail for a mem- 
brane with parameters commonly employed in micro- 
fabrication. 

The silicon nitride optical test films were prepared by 
a LP-CVD (Low-Pressure Chemical- Vapor-Deposition) 
process optimized for low tensile stress and refractive in- 
dex [8]. The 5:1 SiH 2 Cl 2 / NH 3 gas ratio employed results 
in a tensile stress < 100 MPa and optical index greater 
than ~ 2 [9]. The test structure is shown schematically 
in Fig. 1 (inset). Double-side-polished silicon (75- mm- 
diameter, 500-/im-thick) wafers |10| were used as a me- 
chanically robust handling structure for the SiN x mem- 
branes. A 150-nm thermal oxide was grown on the silicon 
wafers by wet oxidation at 950°C for 31 minutes. This 
layer was subsequently used as an etch stop to protect 
the nitride during definition of the silicon handling wafer 
geometry. A low-stress SiN x layer was then deposited by 
LP-CVD (e.g., deposition parameters for 2-/im film are 
835° C for 9.7 hours with pressure 33 Pa and 12 seem 
NH3, 59 seem S1H2CI2). The wafers were then patterned 
with a resist mask and SiN x /Si02 windows formed by 
deep reactive ion etching which removed all the silicon 
under the window area. The residual thermal oxide was 
removed with HF vapor etch leaving a set of uniform 
SiN x membranes each with a 10-mm diameter aperture 
individually suspended from the silicon handling frame. 

The optical tests were performed on SiN x samples hav- 
ing membrane thicknesses of 0.5 and 2.3 fim with a un- 
certainty of 3%. Fabry-Perot resonators were made by 
stacking multiple samples with silicon standoff frames 
between adjacent samples to explore the long-wavelength 
response of the material in greater detail. The silicon 
standoffs allowed a vent path for evacuation of air be- 



tween the nitride membranes. All optical measurements 
were performed in vacuum with a residual pressure less 
than 100 Pa. 

The samples were characterized with a Bruker 125 
high-resolution Fourier Transform Spectrometer (FTS) 
and were measured in transmission at the focal plane of 
an //6 beam. A number of different sources, beam split- 
ters, and detector configurations were used in combina- 
tion to provide measurements over the reported spectral 
range. The single-layer SiN x sample transmission was 
measured over an extended range from 15 to 10000 cm -1 . 
The mercury lamp and a multilayer Mylar beam splitter 
were used to access frequencies below 600 cm -1 . Addi- 
tional mid- infrared spectral data up to 2400 cm~' were 
acquired using a ceramic glow bar source, Ge-coated KBr 
beam splitter, and room-temperature DTGS detector. 
The remaining near-infrared data up to 10000 cm -1 were 
taken with a W filament source, Si on CaF 2 beam split- 
ter, and a liquid-nitrogen-cooled InSb detector (Fig. [TJ. 
Far-infrared data between 15 and 95 cm -1 were taken 
using a mercury arc lamp source and a liquid-helium- 
cooled 4.2-K bolometer. Mylar beam splitters of 50-, 75- 
and 125-^tm thicknesses and a multilayer Mylar beam 
splitter were used during separate scans (Fig. [2|. The 
resultant transmission data were merged into a single 
spectra using a signal-to-noise weighting for subsequent 
parameter extraction. 

The dielectric response is represented as a function of 
frequency, u>, by the classical Maxwell-Hclmholtz-Drude 
dispersion model 
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where M is the number of oscillators and e r = e' r + ie" is 
a complex function of (5M+2) degrees of freedom, which 
are as follows: the contribution to the relative permittiv- 
ity £00 = £m+i of higher lying transitions, the difference 
in relative complex dielectric constant between adjacent 

oscillators Aej = ej — £j+x which serves as a measure of 
the oscillator strength, the oscillator resonance frequency 
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u) T . , and the effective Lorentzian damping coefficient T'j, 
for j = l, M. The following functional form is used to 
specify the damping: 
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where ay allows interpolation between Lorentzian (ay = 
0) and Gaussian wings (ay > 0) similar to the approach 
in [12]. The form indicated above enables a more ac- 
curate representation of relatively strong oscillator fea- 
tures. 
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Fig. 1. (Color online) Room-temperature transmission of 
a silicon nitride sample 0.5 /im thick: measured (grey), 
model (black dotted), and residual (red). The shaded 
band's width delimits the estimated 3cr measurement un- 
certainty. A 30 GHz (1 cm -1 ) resolution is employed for 
the measurement. The insert depicts the geometry of the 
SiN x membrane and micro-machined silicon frame. 

The impedance contrast between free space and the 
thin-film sample forms a Fabry-Perot resonator. The ob- 
served transmission can be modeled |13| as a function of 
the dielectric response (Eq. [T]), thickness, and wavenum- 
ber. The dielectric parameters were solved by means of a 
non-linear least-squares fit of the transmission equation 
to the laboratory FTS data. Specifically, a sequential 
quadratic programming (SQP) method with computa- 
tion of the Jacobian and Hessian matrices 14 15 was 
implemented. The merit function, x 2 ,was used in a con- 
strained minimization over frequency as follows: 
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where N is the number of data points, T the mod- 
eled transmittance, T FTS the measured transmittance 
data, and h is the measured sample thickness. We are 
guided by the Kramers-Kronig relations in defining con- 
straints for a passive material: \ej\ > |ej+i|, £y > 



and e r (0) = £\ 16 . For accurate parameter determi- 



adequately transparent to achieve high signal-to-noise, 
and have diffuse scattering as a sub-dominate process. 
The method requires an a posteriori numerical verifica- 
tion for Kramers-Kronig consistency. In the example pre- 
sented here, a numerical Hilbert transform | 17| of 
reproduces e' r (uj) to within 2% (Fig. [3]). An alternative 
method employing reflectivity and phase allows a priori 
Kramers-Kronig consistent results [18] . However, given 
the details of the thin-film samples and available instru- 
mentation, this approach was not implemented. 
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Fig. 2. (Color online) Measured (solid grey) and model 
(black dotted) transmission for a 3-layer stack of silicon 
nitride samples 2.3 /jm in thickness with 998- /jm inter- 
membrane delays which complements the data shown in 
Fig. [TJ The sample response in the far-infrared was ac- 
quired with a resolution of 3 GHz (0.1 cm -1 ). 

Figure [TJ illustrates the measured and modeled results 
obtained from the analysis of a 0.5-/^m-thick sample. The 
peak residual in the transmittance is less than 3% and 
the 3rj = 0.023 uncertainty band indicated corresponds 
to the 99.7% confidence level. The standard deviation 
adopted for the measured data, a, was estimated assum- 
ing the errors as a function of frequency are uniform 
and have a reduced \ 2 equal to unity. An additional 
uncertainty in the FTS normalization influences the di- 
electric response function at the 1% level. In addition to 
the channel spectra, the observed spectrum shows two 
predominant features at 12 THz and 25 THz. Simula- 
tions with M = 2 oscillators lead to a peak residual 
on transmission of 5% and do not enable recovery of the 
resonance at 25 THz. Using 5 oscillators satisfactorily re- 
covers the observed transmittance and reduces the peak 
residual by a factor of 4.4. When the resonator's quality 
factor, Qeff. — Uj/T'j, is greater than 5, the data were 
not reproducible by either a pure Lorentzian oscillator or 
Eq. (4.6) in [12]. In these regions, the peak transmission 
residuals were decreased by a factor ~ 2 through the use 
of Eq. (2j. 

In Fig. [3] the values of the real and imaginary compo- 
nents of the dielectric function are illustrated as a func- 
tion of frequency. The uncertainty in e r was propagated 
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and computed as described in [l9j . Table [I] contains a 
summary of the best fit parameters for 5 oscillators, 
which can be used to reproduce the data shown in Fig. [3} 

Table 1. Fit parameter summary 



Frequency [cm -1 ] 

10 2 10 3 



10 4 



3 


4 


e" 

3 


uj t . /2ir 


IW27T 




H 




H 


[THz] 


[THz] 


[-) 


1 


7.582 





13.913 


5.810 


0.0001 


2 


6.754 


0.3759 


15.053 


6.436 


0.3427 


3 


6.601 


0.0041 


24.521 


2.751 


0.0006 


4 


5.430 


0.1179 


26.440 


3.482 


0.0002 


5 


4.601 


0.2073 


31.724 


5.948 


0.0080 


6 


4.562 


0.0124 









In order to characterize the long-wavelength portion 
of the dielectric function, Fabry-Perot resonators were 
realized from 1-, 2-, and 3-layer samples. Representa- 
tive data for the 3-layer resonator stack is presented 
in Fig. [2j A multilayer transfer matrix analysis 13 



is 

used to extract the dielectric function using the meas- 
ured SiN^ (2.3 /im) and silicon spacer (998 /zm) thick- 
nesses. The circular symbols at 1.5 THz and 2.5 THz in- 
dicated in Fig. [3] were computed from a composite analy- 
sis of the 3 Fabry-Perot measurement sets. The horizon- 
tal range indicates the data used in each fit. The best 

estimates are e r ~ 7.6 + i0.08 over the range 2-3 THz 

and e r fa 7.6 + iO.OA over 0.4-2 THz. The real compo- 
nent of the static dielectric function derived from the 
data is in agreement with prior reported parameters for 
this stoichiometry [3] . As shown in Fig. [3j the measure- 
ments are internally consistent and represent roughly a 
factor-of-three reduction in uncertainty relative to prior 
infrared SiN x measurements identified by the authors 
[5}j7] . The dielectric parameters reported here are repre- 
sentative of low-stress SiN x membranes encountered in 
our fabrication and test efforts. 
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